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We present a detailed, realistic analysis of the implementation of a proposal for a quantum phase 
gate based on atomic vibrational states, specializing it to neutral rubidium atoms on atom chips. 
We show how to create a double-well potential with static currents on the atom chips, using for 
all relevant parameters values that are achieved with present technology. The potential barrier 
between the two wells can be modified by varying the currents in order to realize a quantum phase 
gate for qubit states encoded in the atomic external degree of freedom. The gate performance is 
analyzed through numerical simulations; the operation time is ~ 10 ms with a performance fidelity 
above 99.9%. For storage of the state between the operations the qubit state can be transferred 
efficiently via Raman transitions to two hyperfine states, where its decoherence is strongly inhibited. 
In addition we discuss the limits imposed by the proximity of the surface to the gate fidelity. 



I. INTRODUCTION 

The idea of encoding information in quantum two-level 
systems (qubits) instead of in classical bits promises a 
revolution in the way we process and communicate in- 
formation 0. Quantum computers, i.e., processing units 
manipulating qubits instead of classical bits, would lead 
to an intrinsic speed-up of calculation that is not possible 
with a classical computer For this purpose, a set 

of controlled operations is necessary, that substitute the 
network of electronic logic gates of present microelectron- 
ics. Analog to classical bits, also for qubits a limited set 
of universal gates exists that allows to implement net- 
works for the execution of any quantum algorithm 0. 
Such universal gates operate on single qubits and on pairs 
of qubits. Whereas the design of single-qubit operations 
is conceptually simple and experimental implementations 
are within the reach of current technology in many cases, 
two-qubit operations arc more demanding in terms of 
both theoretical and experimental investigations. 

Several schemes for the realization of two-qubit gates 
with atoms or ions have been proposed in the last 
years |E S B S ^ £^nd a considerable amount of 
progress has been recently realized for the future imple- 
mentation of quantum information with cold atoms in 
optical lattices [ill IT^ . Many of these schemes have 
been elaborated for ideal systems (harmonic oscillators, 
two-level systems, etc.) or under simplifying approxima- 
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tions for the sake of simplicity, in order to illustrate the 
idea on which the gate is based. The realization of such 
schemes in an experimental setup, where the characteris- 
tics of the real physical system are inevitably to be taken 
into account, is often a formidable task. The inclusion 
of realistic features and values requires a reexamination 
and modification of the schemes presented in the litera- 
ture, in order to exploit all properties of the real physical 
system. For instance, the collisional gate with switching 
potentials proposed in exploits the complete revivals 
of wave packets in harmonic traps, which arc difficult 
to achieve for neutral atoms, and deviations from har- 
monicity pose limiting restrictions to the performance of 
the gate [lj| . More realistic proposals for quantum gates 
are thus required. 

With the spirit just outlined, we reconsider the two- 
qubit gate initially proposed in Ref. |lOl | . Two atoms are 
trapped in a double well potential. The ground (|g)) and 
first excited (|e)) vibrational states of each well are used 
as qubit logic states and the phase gate operation occurs 
via selective interaction between the atoms in the clas- 
sically forbidden region under the potential barrier that 
separates the two atoms. In this paper we investigate 
the feasibility of such a scheme with neutral atoms on 
atom chips. Instead of assuming ideal trapping potentials 
for the atoms to facilitate numerical analysis, we derive 
here the exact potential created by appropriate station- 
ary and time-dependent currents and bias magnetic fields 
on atom chips. The finite size of the current-carrying 
wires is also taken into account and realistic values for 
all physical parameters are used. In addition we consider 
also effects of the proximity of the surface, like spin flips 
and decoherence. Moreover, we specialize our discussion 
to neutral ^^Rb atoms. In this way our investigations 
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are much closer to the experimental conditions and our 
results show how the implementation of such a scheme is 
indeed feasible. 

In spite of the improved realism of our investigations, 
some features of the real system arc neglected in our anal- 
ysis. Nonetheless, we expect that their role is marginal 
and/or can be minimized. Van deer Waals and Casimir 
forces between the atoms and the chip surface are not 
taken into account, since these effects become negligible 
when the atoms arc kept sufficiently distant from the chip 
surface. Similarly, fabrication defects, such as roughness 
of the chip surface and imperfections in the atom chip 
are not taken into account, since they were shown to be 
sufficiently small, and the potentials sufficiently smooth 
when using the appropriate fabrication techniques 
[l5j . Thermal fluctuations of currents in the metal layers 
of the atom chip are another important issue [Til IT^ 
They may lead to loss of the atoms, and to decoherence 
of quantum states. However, we shall show that one can 
transfer the qubit state from the external degree of free- 
dom to the internal one, using two clock states whose 
decoherence is strongly reduced |l9t \2(1^ and the impor- 
tant aspect to consider is the atom loss. Therefore, in 
spite of these approximations, our analysis is a signif- 
icant step towards the first implementation of a phase 
gate with neutral atoms on atom chips. 

The outline of the paper is as follows: in Sec. |n] we 
discuss the properties of ®^Rb atoms that are relevant 
for our analysis and describe the architecture of wires 
and bias magnetic fields on atom chips that create a 
double-well potential for the two atoms with strong con- 
finement in one dimension. We also give the values of 
currents, bias fields, and wire sizes that realize this trap. 
In Sec. mil we give the results of our numerical investi- 
gations on the performance of the phase gate. We show 
that high fidelity 99.9%) can be achieved with short 
operation times 10 ms). In Sec. IIVI we show how to 
transfer the qubit state from the external degree of free- 
dom to two hyperfine (magnetic tappable clock) states 
with Raman transitions. The relative phase of the two 
hyperfine states is insensitive, to first order, to fluctua- 
tions of the magnetic fields, thus, preserving the purity 
of the qubit state for long time. Finally in Sec.|3we re- 
late the achievable gate operation times to the expected 
atomic life times and coherence time, to estimate a real- 
istic fidelity. Whereas the decoherence can be avoided by 
storing the qubits in the clock states, the loss of qubits 
due to thermally induced spin fiips remains, and limits 
the fidelity of operations. The Appendix contains details 
concerning the magnetic field created by the currents on 
the atom chip. 



II. DOUBLE- WELL POTENTIALS WITH 
MAGNETIC FIELDS ON ATOM CHIPS 

Atom chips are versatile integrated microstructures for 
the manipulation of samples of atoms in the ultracold 
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FIG. 1: (Color online) Schematic view of the atom chip config- 
uration. The two wires along the j/-axis lie on the chip surface 
and are separated from the quadrupole wire by 200 nm. 

and quantum degenerate regime (sec Ref. pH and ref- 
erences therein). Trapping potentials for neutral atoms 
can be created near the chip surface. We shall consider 
the simple case of homogeneous bias magnetic fields and 
magnetic fields created by dc (but time-dependent) cur- 
rents. The spin of slow, cold atoms remains adiabatically 
aligned with the magnetic field. The trapping magnetic 
potential, in the weak field approximation, is expressed 

by 

V{r) = -gpUBmpBir) (1) 

where is the Bohr magneton, gp is the Lande factor, 
mp is the azimuthal quantum number, and B(r) is the 
magnetic field. In section IIIII a higher order accurate 
formula - derived from the Brcit-Rabi approach - 
will be used in the numerical simulations, but we limit 
our discussion to the first order term given in Eq. ^ for 
the sake of simplicity. 

The most simple configuration of wires and bias mag- 
netic fields that creates a double-well potential is shown 
in Fig. n] [23. A longitudinal wire along x (hereafter, 
quadrupole wire) carrying a dc current Iq and a uniform 
bias magnetic field B^y perpendicular to the wire create 
a quadrupole potential, with a zero magnetic field along 
a line parallel to the quadrupole wire. Clearly, along this 
line the magnetic field is minimum, but Maiorana spin 
flips occur at zero magnetic field, which result in atom 
losses. Therefore the minimum is shifted to a non-zero 
value with the addition of a second uniform bias mag- 
netic field Bqx, orthogonal to the first one and parallel 
to the chip surface. Two more wires (hereafter, left and 
right wire, respectively), perpendicular to the quadrupole 
wire, carry a dc current I1.2 — ctlo, whose magnetic fields 
give rise to a modulation of the trapping potential. 

In consideration of the small extension of the system 
(a few microns), the assumption of infinitely thin wires 
is too crude an approximation. Therefore in the poten- 
tial we shall use the expression of the magnetic field 
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created by currents flowing through wires with finite rect- 
angular sections (a typical shape of wires on atom chips). 
Details of the calculation are given in the Appendix. 

A trapping potential with two well-separated min- 
ima, as shown in Figure l^a), is created with the val- 
ues Jo = 40.89 mA, a = 70.25 x lO^^, Bq^ = -9.90 G, 
and Boy = 50.0 G. The three wires have an identical 
rectangular section, of width W = 700 nm and height 
H = 200 nm. and the centers of the left and right wires 
arc 1.60 /^m apart. The center of the quadrupole wire 
is at a distance zq = 400 nm under the chip surface, 
whereas the left and right wires lie on the chip. We stress 
that the values for the required transversal confinement 
has been approached even in the first experiments with 
nano fabricated atom chips [^13 ^'^d the currents, bias 
fields, size and distances of the wires are within current 
laboratory conditions |2^ . 

The two potential minima are at a distance of 1.19 ^m 
from the surface, and the line joining them is slightly 
tilted by an angle (3 ~ 14.8° from the x axis. This an- 
gle defines the new axis x' along which the dynamics 
will take place (see |2^ for details). We also define a 
new axis y' parallel to the chip surface and perpendic- 
ular to x' . The z axis remains unchanged. Since the 
trapping frequencies at the two minima verify the con- 
dition Vx' <C Vyi — Vz p^ . transverse vibrations do not 
get excited during the gate operation (provided that the 
energy of the atoms is less than the quanta of trans- 
verse vibrations), and the trapped atoms will experience 
a quasi one-dimensional (ID) dynamics. 

The value of the magnetic field at the two minima is 
Bm — 3.23 G. This value minimizes the decohcrencc in- 
duced by fiuctuations of the dc currents for the hyper- 
fine states \F = 2,mp = 1) and \F = l,mF = —1) of the 
ground state 5Si/2 of ^^Rb These clock states will 

be used to store the qubit information at the end of the 
gate operation, as described in Sec. IIVI 



III. GATE OPERATION 

Compared to previous static proposals j2^ . we imple- 
ment here a time-dependent phase gate in a spirit similar 
to the one of Ref . [ij • 

As it can be noticed in Figure |2Ia), when the barrier 
is high the translational wavefunctions of the atoms do 
not overlap in the inter-well region. In this type of envi- 
ronment, the atoms do not interact. On the other hand, 
when the barrier is lowered, as in Figure |2Ib). tunneling 
takes place and the probability of finding the atoms in 
the classically forbidden region is not negligible any more. 
As a consequence, the energy splitting between the sym- 
metric and anti-symmetric state combinations increases 
quickly when the barrier height ^ is lowered. This ef- 
fect is clearly state-selective in the sense that it affects 
differently the ground and excited translational states. 
It therefore constitutes an interesting candidate for the 
implementation of a conditional logical gate. 
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FIG. 2: (Color online) Double well potentials created by the 
atom chip configuration shown in Figure Q The energies 
of the first six eigenstates are shown as red horizontal lines. 
The blue dashed line represents the wavefunction of the third 
eigenstate labeled as Is because it originates from the sym- 
metric combination of the v = 1 trapped levels, also labeled 
as |e) in the text, (a) Highest barrier ^/2n = 35.4 kHz ob- 
tained with /() = 40.89 mA and a = 70.25 x 10"^. (b) Lowest 
barrier = 14.4 kHz obtained with Jo = 42.01 mA and 

a = 69.70 X 10"'^. In both cases the bias magnetic fields are 
equal to Bo^ = -9.90 G, and Boy = 50.0 G. 




FIG. 3: (Color online) (a) Variation of the dc cur- 
rent Io{t) (red lines) in the quadrupole wire and of 
a(t) = Ii{t)/Io(t) = l2{t)/Ioit) (blue lines) during the gate 
operation. The solid lines correspond to a simple linear gate 
and the dashed lines represent an optimized gate (see text for 
details), (b) Variation of the barrier height ^(t) (green dashed 
line) and of the energies of the first six instantaneous eigen- 
states of the double well potential (solid lines) in the case of 
the optimized gate. The times To, (To+Ti) and (2To-|-ri) de- 
limit the three steps which constitute the conditional n phase 
gate. 



In the present scheme, the barrier height £^{t) is con- 
trolled by varying simultaneously the intensities Io(t) and 
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Ii{t) = hit) = a(t) X Io{t) in the quadmpole and in the 
perpendicular wires. In a first and simple implementa- 
tion of the phase gate, we impose a linear variation of the 
barrier height ^ with time. The phase gate is decomposed 
in three steps; 

• When ^ t ^ To the barrier is lowered and the 
double-well potential changes from the one of Fig- 
ure |2I a) to the one of Figure |2Ib). 

• When Tq ^Tq + Ti the inter-well barrier is 
fixed at its lowest value, such that a large inter- 
atomic interaction takes place. 

• Finally, when Tq + Ti t 2To + Ti the inter- 
well barrier is raised again until the initial condition 
is recovered. 

The linear variation of ^(t) is obtained by changing Io{t) 
and a{t) only, as depicted by the solid lines shown in 
Figure Ela). We have verified that this simultaneous 
variation of the dc currents does not modify the direc- 
tion x' along which the dynamics is taking place. In 
addition, the value of the magnetic field at the potential 
minima remains equal to 3.23 G during the whole gate 
operation. The dynamics is therefore quasi-adiabatic if 
Tq^ l/vx' — 77 fis (hereafter we write x,y instead of 

The gate operation is followed by solving the time- 
dependent Schrodinger equation along the double-well 
direction x for the wave packet '^{xi,X2,t) describing 
the motion of the two atoms 

d 

ih—^[xi,X2,t) ^n(xi,X2,t)<^{xi,X2,t) . (2) 

The two-dimensional time-dependent Hamiltonian is 
written as 

H{xi,X2,t) = Tj;j +fx2 +V2Bixi,X2,t) , (3) 

where Tq denotes the kinetic energy operator along the 
g-coordinatc. The two-dimensional potential is given by 
the following sum 

V2B{xi,X2,t) = V{xi,t) + V{x2,t) 

+ Vint{\x2~Xi\,t), (4) 

where V{x, t) is the trapping potential (QJ created by the 
atom chip and T4nt(|a^2 — xi\,t) represents the averaged 
interaction potential between the two atoms at time t 

1/2 

Vint{\x2 - Xi\,t) ^ 2hao{uJyUJz) ' S{\x2~Xi\). (5) 

This last expression is obtained by averaging the three- 
dimensional delta function interaction potential over the 
lowest trap states along the y and z directions 0- One 
can note that the atom-atom interaction strength is pro- 
portional to the s-wave scattering length uq. Since the 
orthogonal trapping frequencies tOy and u)z vary slightly 
during the gate operation |2^ . the averaged interaction 
strength is also slightly time-dependent. 




FIG. 4: (Color online) Infidelity of the conditional phase gate 
as a function of the gate duration. The red and blue lines 
correspond to the infidelities calculated for the linear and op- 
timized gates respectively (see Figure |3 and text for details). 

We solve the time-dependent Schrodinger equation (0 
in a basis set approach by propagating the initial state of 
the two-atom system in time. We start with the atoms 
initially in one of the first four eigenstates \gg), Iffe), \eg), 
or |ee) of the double-well potential 10} . The wavefunction 
^(xi,X2,t) is then expanded as 

*(xi,a;2,t) = ^Ci(t) X (pi{xi,X2,t) (6) 

i 

where ipi{xi,X2,t) represent the wavefunctions asso- 
ciated to the instantaneous eigenstates of the two- 
dimensional potential V2D {xi ,X2,t) of Eq. Q . Insert- 
ing this expansion into the time-dependent Schrodinger 
equation ^ yields the following set of first order cou- 
pled ordinary differential equations for the complex coef- 
ficients Ci{t) 

ih^c^it) = £i Ci{t) - ih^Cj{t) Vtj{t) , (7) 

j 

where Si denotes the energy of the eigenstate (pi and 
Vij (t) is a time-dependent non-adiabatic coupling aris- 
ing from the time variation of the barrier height ^(t) 

V^Jit) = {ip,\ - \ip,) = ^ (^,1 — l^j) . (8) 

This set of equations is solved using an accurate 
Shampine-Gordon algorithm 12^. For each value of the 
duration Tq we adapt Ti in order to obtain an accurate 
conditional phase gate, with = 7r±lO~^. 
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At the end of the propagation {t = tf) the coefficients 
Ci{tf) are analyzed to calculate the infidelity of the gate 



/ = 



E 



--\gg)- 



^itf)\■ 



(9) 



where the sum runs over all possible initial qubit states. 
The infidelity is therefore a measure of the deviation from 
adiabaticity which arises from the non-adiabatic cou- 
plings Vij (t). This quantity is plotted in Figure^ (red hue 
for the linear gate) as a function of the gate duration. It 
shows an oscillatory behavior partially similar to the one 
observed with atoms trapped in an optical lattice |lCll| . 
The succession of maxima and minima is a signature of 
constructive and destructive interferences between two 
distinct pathways of excitation of the initial qubit state. 
Indeed the initial state may be excited in the time in- 
tervals ^ t ^ To and To + Ti t 2To + Ti, when the 
barrier is lowered and raised. The nature of this inter- 
ference depends on the phases which develop during the 
gate operation The periodicity of the oscillation is 
simply related to the Bohr frequencies associated with 
the energy splitting of the two-atom eigenstates. The 
linear gate configuration proposed here can achieve a rel- 
atively high fidelity of about 99.6% in just 7.6 ms. 

One should also realize that in the general case the 
couplings between the initial qubit states and the other 
accessible two-atom eigenstates vary with time. These 
couplings effectively increase when the inter-well barrier 
approaches the energy of the initial state. The linear gate 
proposed until now is therefore far from being optimal for 
the maximization of the gate fidelity. 

We have thus implemented an optimized gate which 
tends to minimize these couplings during the whole gate 
duration. For this purpose, we impose a fast variation 
of the barrier height £,{t) at early times i <C Tq, whereas 
this variation is much slower when i ~ Tg- This is done 
by choosing 



at 



Si - e 



\ee) 



(10) 



In this expression, 7 is a dimensionless proportionality 
factor, which can be decreased to achieve larger gate 
durations. The first derivative with respect to time of 
the barrier height ^(<) is therefore chosen such that the 
maximum effective first-order coupling between the high- 
est energy qubit state |ee) and all other states remains 
constant during the whole gate duration. With this ap- 
proach, higher fidelities are expected when compared to 
a linear gate of the same duration. 

The variations of Io{t) and a{t) for this optimized gate 
are shown as dashed lines in Figure |2Ia). Figure 0] shows 
that the optimized gate infidelity (blue line) is, on aver- 
age, improved by a factor 6 when compared to the linear 
gate. As a consequence, this optimized gate can achieve 
fidelities of 99% in 6.3 ms and of 99.9% in 10.3 ms. 



IV. USING BOTH DEGREES OF FREEDOM AS 
QUBIT STATES 

When neutral atoms are used, it is possible to em- 
ploy either the vibrational or the internal states as qubit 
states. Section Hill has shown that vibrational states are 
very promising in terms of gate performance. However, 
the readout of the qubit after the gate operation seems to 
be a difficult task in comparison to the readout of inter- 
nal states, for which observation of fluorescence emitted 
from selected transitions is a well established technique. 
Moreover, a closer examination of the hyperfine structure 
of the ground state of ^^Rb provides further motivation 
for the use of internal states. 

The hyperfine structure of the 5Si/2 ground state 
of ^^Rb is shown in Figure O eight sublevels have 
to be considered, three with = 1 and five 

with F = 2. Only three of these eight sublevels 
(the states |F = 1, mij- = — 1), \F = 2,mp = l), and 
\F =■ 2,mp = 2)) are low-field seeking states which can 
be trapped by static magnetic potentials [23|. In ad- 
dition, the two hyperfine levels \F = X^mp = —\) and 
\F = 2,mp = 1) have opposite Landc factors. As a con- 
sequence, they experience identical magnetic potentials. 
A more precise estimation that goes beyond the linear 
approximation given in Eq. fQ) shows that the difference 
of the Zeeman shifts experienced by these two states is 
minimum aX B ~ 3.23 G At this field the states 

|0) = \F = 2,mF = 1) and |1) = \F = l,mp = -1) are 
insensitive to field fluctuations at first order. The deco- 
herence of an arbitrary superposition of the two states 
|0) and |1) due to fluctuations in the current intensities 
(which entail fluctuations in the magnetic field and thus 
in the potential CQl) is thus strongly inhibited. Indeed, 
coherent oscillations between the states |0) and |1) have 
recently been observed and the decoherence time has 
been estimated to be Tc = 2.8 s ±1.6 s The inhi- 

bition of decoherence is a fundamental issue in quantum 
computation, since the coherence of the quantum state is 
an essential ingredient. It would thus be desirable to use 
the clock states \F = l,mp = —1) and \F = 2,mp = 1) 
for the storage of quantum information. 

In fact, we can use both degrees of freedom as qubits: 
when information must be processed, it will be encoded 
in the external degree of freedom and the gate operation 
will be done as described in the previous Section. At 
the end of the gate operation, the qubit state can be 
transferred to the internal degree of freedom and stored 
there until new operations are performed or until the 
final readout. In order to realize this transfer from one 
degree of freedom to the other, we need to swap the states 
associated with these two degrees of freedom, according 
to 



I.9I) ^ |e0) 



(11) 



In this way, an arbitrary qubit state encoded in the vi- 
brational states as 



{a\99) + b\ge) + c\eg) + d\ee)) (g, |00) 



(12) 
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FIG. 5: (Color online) Schematic representation of the hyper- 
fine structure of the 5 S1/2 and 5 P1/2 states of *^Rb . The two 
qubits |0) = \F = 2,m.F = 1) and |1) = \F = l,mF = -1) of 
the 5Si/2 ground state are shown in red. The intermediate 
sublevels |A), |B), |C), and |D) of the 5Pi/2 excited state are 
shown in blue. The cr+ and g- transitions connecting them 
are shown with dotted (orange) and dash-dotted (green) ar- 
rows respectively. 



is transformed into 

\gg) ® (a|00) + &|01) + c|10) + d|ll)) (13) 

and vice versa. Eq. (|12|) assumes that the internal degree 
of freedom of the two atoms has been initialized to the 
state |0). The map given in Eq. (|11|) is a selective switch 
from the internal state |0) ( respectively |1)) to the state 
|1) ( respectively |0)) under the control of the vibrational 
state. In order to realize it, we propose the use of two- 
photon Raman transitions. 

The realization of such a two-photon Raman transition 
via intermediate 5 P states requires a careful analysis (2^ . 
The 5Si/2 and 5Pi/2 levels have eight hyperfine states 
with F = 1 or 2, whereas the 5 P3/2 level comprise sixteen 
states with F = 0, 1, 2 or 3. The Raman transition re- 
quires two lasers, one with polarization (frequency uj+ 
and wave vector k-|_), the other with a~ polarization (fre- 
quency w_ and wave vector k_ ). When the laser frequen- 
cies w+ and u)- are kept below the 5 S1/2 — '■S P1/2 transi- 
tion, the states in 5 P3/2 can be neglected, being detuned 
by about 7 THz. From now on we therefore restrict our 
analysis to the states |0) and |1), and to the following 
5Pi/2 states 
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Taking into account all possible transitions according 
to the selection rules imposed by the polarizations, the 



Hamiltonian in the usual rotating-wave approximation is 



with 
Ho 



H. = Tio + "Hi + Ti.2 



D 

h J2 ^k\k){k\+HY,^j\j){j\ 

k=OS j=A 

ha 



(14) 



H2 = 



2V3 



''+(|A)(1| + |C)(1|-a/2|D)(0| 
(^/2|B)(1| + |C)(0|-|A)(0| 



h.c. 
h.c. 



where we have defined the phases 



V- 



k+ • r 
k • r 



■ CiJ+ t 
LU^ t 



and where ri+ and f2_ are the Rabi frequencies asso- 
ciated with the interaction between the electric dipole 
of the atom and the electric field of the lasers, r is the 
atomic position operator, and h.c. denotes the Hermitian 
conjugate. 

From the Hamiltonian H14|) one can derive the Heisen- 



berg equations of motion for the operators cr^ 



^XjI, 



with z, J = 0, 1, A, 
(TAi we get 



D. For instance, for the operator 



i h&Ai 



[dAi , Ha] + [<7Ai , ni] + [aAi , (15) 



where 



[o-Al, Wo] 
[o'Al,'Wl] 

[o'Al,'W2] 



h {uji - lja) 0-Ai 



2V3 



CTAA - CTii + (TAG , 



V2 



CAB 



With the introduction of the "slow" variables 

i Lj^ t 



o^Ai = CTAi e 
o'oi = ctqi e 
equation H15|) becomes 



-i —LJ^ ) t 



I (TAI 



{UJ+ -UJX+UJl) CTAl 



2V3 



^ e I (TAA - CTll 



CTAC 



) 



I O'OI + V2e ^ +' (TAB I 



When the lasers are strongly detuned from the atomic 
transition frequencies, one can obtain an effective 
Hamiltonian by imposing the adiabaticity condition 
(TAi — [33,133. In this case, we get 



CTAl 



+ e 



2^3 Aai 

Vt^ e-"'i- 



C^AA - 0-11 ^- (Tac^ 

+ V2(tabV 
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where we have defined the detuning 



(16) 



Performing similar calculations for the operators ctbi, 
ccij CTAo, Ceo and (Tdo after substituting them into Hi 
and 7^2 and including the diagonal terms into Ho, we get 
the effective interaction Hamiltonians 



Hi 



2V6 



17 _ e "' 



C^AB 



/ ^AC 



O'AC 



V2 



Aai/ ' ^/6 V^Ai ' Aci 

CAD , CTbc CCD 



Ado 
coi 



Aci 

coi 

Aci 



Ado 
+ h.c 



and 

H2 



(17) 



+ ^ e 



CAB \ _ / CAC CAC 

A^y 72 V aI^^ a^ 

, Cbc CCD 



Cad 

Aao 
coi 



Abi 
coi 



Aco 
f h.c. 



(18) 



a/6 V Aao Aco , 
where the phase 77 is simply given by the difference 

V = V+- V- (19) 
and the detunings are defined just as (|16|l by 



Aao = ujA - ujQ - uj^ 
Aco = wc - ^0 - tc)- 
Ado ^ t^D — ^0 ^ ^+ 



Aai = WA - - 

Aci = cjc — t^i - 
Abi = cjb - ti'i - ^- 



When the initial atomic state is a superposition of |0) 
and |1), as in our case, we can retain the terms contain- 
ing (jQi and CTio only, such that the effective interaction 
Hamiltonian reduces to 



H'"' = H[ + H'2 



where 



4V3 



Cqi 

aI7 

Coi 

aI^ 



Cqi 
A^ 
Cqi 
A^ 



(20) 

+ h.c. (21a) 
+ h.c. (21b) 



The presence of two 5Pi/2 states with rap = (the 
I A) and |C) states) could seem problematic. When the 
lasers are very far detuned from resonance, the detun- 
ings from these two levels are almost equal: Aao — Aco 
and Aai — Aci. A destructive interference between 
these two quantum paths suppresses the Raman tran- 
sition between the states |0) and |1) (see Eqs. (j2D)- 
In order to stimulate this Raman transition, the detun- 
ings from levels |A) and |C) must be significantly differ- 
ent. This condition can be realized in ®^Rb , where the 
states |A) = \F = l,mF = 0) and |C) = = 2, tof = 0) 



are separated by loq — lja — 27r • 800 MHz and the natu- 
ral line width of the level 5Pi/2 is about 27r • 5.75 MHz. 
Finally, we comment the realization of the sideband exci- 
tation required by the transformation (|ll|l . It requires a 
laser linewidth much smaller than the vibrational angular 
frequency, i.e., much smaller than 10 KHz, a condition 
that can be reached experimentally with standard tech- 
niques 01 . It is also possible to produce pulses with spe- 
cific shapes that optimize the transition probability and 
stimulate the transition in very short periods, so that the 
gate operation time is not significantly altered. 



V. LOSS AND DECOHERENCE 

An important mechanism we have to consider in a re- 
alistic evaluation of the performance of quantum gates, is 
the possibility of loss and decoherence of the qubits dur- 
ing the operation. Especially in atom chips this is an im- 
portant issue, since the cold atoms are perturbed by ther- 
mal electromagnetic fields generated by the nearby, hot 
solid substrate. This leads to heating, trap loss and deco- 
herence 0, 0, 01 . Since the thermal induced couplings 
lead to similar timescales for all these undesired processes 
we consider the loss, which was studied in largest detail 
theoretical and experimentally |33 |. In addition, as we 
have seen in the previous section, the decoherence can be 
dramatically reduced by transferring the qubits o clock 
states during idle periods. 

Following the treatment of one can estimate the 
atomic lifetimes for a setup like shown in Fig.^ Assum- 
ing gold wires we find a typical lifetime of 0.8 s due to 
thermally induced spin flips and consequent loss of the 
qubit state. The left and right wires, which are closest 
to the atoms and carry small currents, give the largest 
contribution to the loss. In this configuration the best 
fidelity of ^ 98.7% is obtained for a gate operation time 
of ~ 9 ms. 

Small changes in the trapping geometry can improve 
the fidelity: Since a big part of the loss comes from the 
left and right wires, one can fabricate them much thinner. 
Reducing the thickness of the double well wires to 50 nm 
would increase the lifetime to over 2 sec and increase 
the fidelity to ~ 99.5 % at a gate operation time of ~ 
10 ms. This can be even more improved by a little bit 
different wire geometry for the quadrupole wire. A wider 
but thinner quadrupole wire (~ 3/im x 50nm), and 
smaller left and right wires 0.3/im x 50nm) will lead 
to lifetimes in excess of 3 s and consequently a better 
gate fidelity in the order of ~ 99.7 % at a gate operation 
time of ^ 11 ms. 

Further improvements can be expected when optimiz- 
ing the wire geometries. Two aspects have to be taken 
into account. 

• The best achievable lifetime increases with the 
square of the magnetic field gradient dB/dZ [33l |. 
The present trap has a transverse gradient of 300 



8 



kG/cm. For 100 kG/cm and optimized wire cross 
section one can achieve lifetimes of up to 100 s |33j . 

• A second consideration is using different materials, 
and cooling the chip surface, as investigated theo- 
retically by 1131 ■ This should lead to even further 
improvements. 

• Directly fabricating the wires into semiconductor 
chips will give even better thermal coupling and 
allow the use of high resistivity materials for the 
wires, which should improve even further the gate 
performance. 

Overall we conclude that for optimized geometries fi- 
delities better than 99.9% should be possible in realistic 
settings with present day atom chip technology. 
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VIII. APPENDIX 



VI. SUMMARY AND CONCLUSIONS 

We have presented a detailed analysis of the imple- 
mentation of a quantum phase gate with neutral rubid- 
ium atoms on atom chips. Our analysis is quite close to 
the experimental conditions and is within the reach of 
current technology. We have shown how to create a dou- 
ble well potential near the surface of an atom chip and 
studied the performance of a phase gate, using as qubit 
states the vibrational states of two rubidium atoms, each 
sitting in one of the two wells. The gate operation is re- 
alized through an adiabatic modification of the potential 
barrier that separates the two wells. This implementa- 
tion allows selective interaction between the vibrational 
states of the two atoms. We have found that, neglect- 
ing losses due to spin flips, a fidehty of 99.9% can be 
achieved in just 10.3 ms. We have also shown that the 
motional qubit state can be transferred to two hyperfine 
states with reduced dccohercnce. The estimation of loss 
due to thermally induced spin flips reduces slightly the 
gate fidelity, and improvements to reduce these effects 
have been discussed. 

The results presented here are a significative improve- 
ment when compared to an implementation using a static 
trap [2^ . For a fidelity above 99% the operation time is 
diminished from about 16.2 ms to 6.3 ms. The opera- 
tion time could be further reduced using non-adiabatic 
changes of the barrier, as recent studies employing opti- 
mal control theory indicate [s^ , but the excitation of the 
perpendicular degrees of freedom might then become an 
issue. 
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The magnetic field created by an infinitely thin and 
infinitely long straight wire is 



B = 



kol 



(22) 



where / is the dc current flowing through the wire, r 
is the distance to the wire and feg = ^o/27r {fJ-o is the 
magnetic permeability of the vacuum). In the case of a 
wire of finite rectangular section, we average the mag- 
netic fields created by many infinitely thin wires lying 
inside the finite section. For instance, for a dc current 
/ flowing in the positive x direction through a wire with 
section {~-W/2 < y < W/2, -H/2 < z < H/2), the 
magnetic field components are 



Byiy^^) = TF77? I dy' I dz' 



HWj_w% J_H/2 {y - y'? + z') 



Ta (23) 



and 



B.iy.z) = T^I dy' I dz' 



y_l (24) 



(y - y'? + {z- z') 

In Eq. H23fl the integration over y' gives the y-component 



Byyy^^i^j^i dz [&Yct&ii- 



y- 



-H/2 



-arctan— ^i— 1 (25) 



where y± = y ^ W/2. A final integration over gives 
kol 



By{y,z) = - 



HW 



z- I arctan arctan — 

y^ y- 



z+ z+ 
arctan arctan — 

y- y+ 

y- , y- + y+ . yl + zl, 
— log , 3 - — log 
2 y-L+z% 2 y+ + z^ 



(26) 
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After defining z± = z zk H/2, analogous calculations for 
the component Bz{y,z) give 



kol \ ( ^ y- . y- 



HW 

+y+ 



?/_ arctan arctan - 

i arctan — — arctan — 

z_ +yl z+ zl+y^_ 

' 2 zl+yl 2 zl + yl 



(27) 
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